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Abstract 

This is the extensive follow-up report of a recent Letter in which the existence 
of self-organized criticality (SOC) in systems of interacting soft gluons is proposed, 
and its consequences for inelastic diffractive scattering processes are discussed. It 
is pointed out, that color-singlet gluon-clusters can be formed in hadrons as a con- 
sequence of SOC in systems of interacting soft gluons, and that the properties of 
such spatiotemporal complexities can be probed experimentally by examing inelastic 
diffractive scattering. Theoretical arguments and experimental evidences support- 
ing the proposed picture are presented — together with the result of a systematic 
analysis of the existing data for inelastic diffractive scattering processes performed at 
different incident energies, and/or by using different beam-particles. It is shown in 
particular that the size- and the lifetime-distributions of such gluon-clusters can be 
directly extracted from the data, and the obtained results exhibit universal power- 
law behaviors — in accordance with the expected SOC-fingerprints. As further con- 
sequences of SOC in systems of interacting soft gluons, the t-dependence and the 
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(M^/s)-dependence of the double differential cross-sections for inelastic diffractive 
scattering off proton-target are discussed. Here t stands for the four-momentum- 
transfer squared, M x for the missing mass, and \/s for the total c.m.s. energy. It 
is shown, that the space-time properties of the color-singlet gluon-clusters due to 
SOC, discussed above, lead to simple analytical formulae for d 2 a / dt d(M% / s) and for 
da/dt, and that the obtained results are in good agreement with the existing data. 
Further experiments are suggested. 



1. Interacting soft gluons in the small-x^ region of DIS 

A number of striking phenomena have been observed in recent deep-inelastic electron-proton 
scattering (DIS) experiments in the small-x^ region. In particular it is seen, that the 
contribution of the gluons dominates!], and that large-rapidity-gap (LRG) events existEri. 
The latter shows that the virtual photons in such processes may encounter "colorless objects" 
originating from the proton. 

The existence of LRG events in these and otherBi scattering processes have attracted 
much attention, and there has been much discussionEHO] on problems associated with the 
origin and/or the properties of such "colorless objects". Reactions in which "exchange" of 
such "colorless objects" dominate are known in the literatureiSi'i as "diffractive scattering 
processes". While the concepts and methods used by different authors in describing such 
processes are in general very much different from one another, all the authors (experimen- 
talists as well as theorists) seem to agree on the following! (see also Refs. ||-f§ ^0)-Q): 
(a) Interacting soft gluons play a dominating role in understanding the phenomena in the 
small-xs region of DIS in general, and in describing the properties of LRG events in par- 
ticular, (b) Perturbative QCD should be, and can be, used to describe the LRG events 
associated with high transverse-momentum (p±) jets which have been observed at HERA0 
and at the Tevatronl Such events are, however, rather rare. For the description of the bulk 
of LRG events, concepts and methods beyond the perturbative QCD, for example Pomeron 
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Models! based on Regge Phenomenology, are needed. It has been suggested a long time ago 
(see the first two papers in Ref. ||) that, in the QCD language, "Pomeron-exchange" can be 
interpreted as "exchange of two or more gluons" and that such results can be obtained by 
calculating the corresponding Feynman diagrams. It is generally felt that non-perturbative 
methods should be useful in understanding "the small-x^ phenomena", but the question, 
whether or how perturbative QCD plays a role in such non-perturbative approaches does 
not have an unique answer. 

In a recent LetterB, we proposed that the "colorless objects" which play the dominating 
role in LRG events are color-singlet gluon-clusters due to self-organized criticality, and that 
optical-geometrical concepts and methods are useful in examing the space-time properties 
of such objects. 

The proposed pictured! is based on the following observation: In a system of soft gluons 
whose interactions are not negligible, gluons can be emitted and/or absorbed at any time 
and everywhere in the system due to color-interactions between the members of the system 
as well as due to color- interactions of the members with gluons and/or quarks and antiquarks 
outside the system. In this connection it is important to keep in mind that gluons interact 
directly with gluons and that the number of gluons in a system is not a conserved quantity. 
Furthermore, since in systems of interacting soft-gluons the "running-coupling-constant" is 
in general greater than unity, non-perturbative methods are needed to describe the local 
interactions associated with such systems. That is, such systems are in general extremely 
complicated, they are not only too complicated (at least for us) to take the details of local 
interactions into account (for example by describing the reaction mechanisms in terms of 
Feynman diagrams), but also too complicated to apply well-known concepts and methods 
in conventional Equilibrium Statistical Mechanics. In fact, the accumulated empirical facts 
about LRG events and the basic properties of gluons prescribed by the QCD are forcing us 
to accept the following picture for such systems: 

A system of interacting soft gluons can be, and should be considered as an open, dynam- 
ical, complex system with many degrees of freedom, which is in general far from equilibrium. 



In our search for an appropriate method to deal with such complex systems, we are led 
to the following questions: Do we see comparable complex systems in Nature? If yes, what 
are the characteristic features of such systems, and what can we learn by studying such 
systems? 

2. Characteristic features of open dynamical complex systems 

Open, dynamical, complex systems which are in general far from equilibrium are not difficult 
to find in Nature — at least not in the macroscopic world! Such systems have been studied, 
and in particular the following have been observed by Bak, Tang and Wiesenfeld (BTW) 
some time ag This kind of complex systems may evolve to self-organized critical states 
which lead to fluctuations extending over all length- and time-scales, and that such fluc- 
tuations manifest themselves in form of spatial and temporal power-law scaling behaviors 
showing properties associated with fractal structure and flicker noise, respectively. 

To be more precise, BTW0 and many other authors^ proposed, and demonstrated by 
numerical simulations, the following: Open, dynamical, complex systems of locally interact- 
ing objects which are in general far from equilibrium can evolve into self-organized structures 
of states which are barely stable. A local perturbation of a critical state may "propagate", 
in the sense that it spreads to (some) nearest neighbors, and than to the next-nearest neigh- 
bors, and so on in a "domino effect" over all length scales, the size of such an "avalanche" 
can be as large as the entire system. Such a "domino effect" eventually terminates after 
a total time T, having reached a final amount of dissipative energy and having effected a 
total spatial extension S. The quantity S is called by BTW the "size", and the quantity 
T the "lifetime" of the avalanche — named by BTW a "cluster" (hereafter referred to as 
BTW-cluster or BTW-avalanche). As we shall see in more details later on, it is of consid- 
erable importance to note that a BTW-cluster cannot, and should not be identified with a 
cluster in the usual sense. It is an avalanche, not a static object with a fixed structure which 
remains unchanged until it decays after a time-interval (known as the lifetime in the usual 
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sense) . 

In fact, it has been shownllMl that the distribution (Dg) of the "size" (which is a measure 
of the dissipative energy, 5) and the distribution (Dt) of the lifetime (T) of BTW-clusters 
in such open, dynamical, complex systems obey power-laws: 

D S (S) ~ S-», (1) 

D T (T) ~ T-\ (2) 

where /i and v are positive real constants. Such spatial and temporal power-law scaling 
behaviors can be, and have been, considered as the universal signals — the "fingerprints" - 
of the locally perturbed self-organized critical states in such systems. It is expectedHEl that 
the general concept of self-organized criticality (SOC), which is complementary to chaos, 
may be the underlying concept for temporal and spatial scaling in a wide class of open, 
n on- equilibrium, complex systems — although it is not yet known how the exponents of 
such power laws can be calculated analytically from fundamental theories such that for 
gravitation or that for electromagnetism. 

SOC has been observed experimentally in a large number of open, dynamical, complex 
systems in non-equilibrium@Hiil among which the following examples are of particular inter- 
est, because they illuminate several aspects of SOC which are relevant for the discussion in 
this paper. 

First, the well-known Gutenberg- Richter lawffl for earthquakes special case of 
Eq. ([[]): In this case, earthquakes are BTW-clusters due to SOC. Here, S stands for the 
released energy (the magnitude) of the observed earthquakes. D S (S) is the number of 
earthquakes at which an energy S is released. Such a simple law is known to be valid for all 
earthquakes, large (up to 8 or 9 in Richter scale) or small! We note, the power-law behavior 
given by the Gutenberg-Richter law implies in particular the following. The question "How 
large is a typical earthquake?" does not make sense! 

Second, the sandpile experiment s00 which show the simple regularities mentioned in 
Eqs.([l]) and @: In this example, we see how local perturbation can be caused by the 



addition of one grain of sand (note that we are dealing with an open system!). Here, we can 
also see how the propagation of perturbation in form of "domino effect" takes place, and 
develops into BTW-clusters/avalanches of all possible sizes and durations. The size- and 
duration-distributions are given by Eqs.(^J) and (0), respectively. This example is indeed a 
very attractive one, not only because such experiments can be, and have been performed in 
laboratories^, but also because they can be readily simulated on a 

Furthermore, it has been pointed out, and demonstrated by simple mo dels0£Hl!, that 
the concept of SOC can also be applied to Biological Sciences. It is amazing to see how 
phenomena as complicated as Life and Evolution can be simulated by simple models such 
as the "Game of Life"0 and the "Evolution Model" Ellll. 

Having seen that systems of interacting soft-gluons are open, dynamical, complex sys- 
tems, and that a wide class of open, dynamical, complex systems in the macroscopic world 
evolve into self-organized critical states which lead to fluctuations extending over all length- 
and time-scales, it seems natural to ask the following: Can such states and such fluctuations 
also exist in the microscopic world — on the level of quarks and gluons? In particular: Can 
SOC be the dynamical origin of color-singlet gluon-clusters which play the dominating role 
in inelastic diffractive scattering processes? 



3. SOC in inelastic diffractive scattering processes? 

One of the main goals of the present paper is to answer the questions mentioned at the end 
of the last section. We discuss in this and in the following four sections, how to look for 
signals of SOC in systems of interacting gluons, and what can we see when we look for signals 
of SOC in scattering processes in which systems of interacting gluons play the dominating 
role. 

Here, we explicitly see: (a) The fundamental properties of the gluons, (b) the necessary 
conditions for the occurrence of SOC, and (c) the available technical possibilities strongly 
suggest that the most favorable place to study the possible existence of such signals is (i) 
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to look at the experimental results obtained in events associated with large rapidity gaps 
in deep-inelastic scattering, (ii) to look at the experimental results in inelastic diffractive 
hadron-hadron scattering, and (iii) to compare such observations with each other. 

Having the special role played by "the colorless objects" in inelastic diffractive scattering 
in mind, let us begin our discussion with the question: What are such "colorless objects"? 
Up to now, we do not know much about such objects. We know that they carry neither 
color- nor any flavor-quantum numbers. We know that they exist in high-energy reactions 
where soft gluons play the dominating role. We know that they can be probed in diffractive 
scattering processes, in the sense that they can interact with different beam-particles. But, 
there are a lot more which we do not know, for example: What is the mass of a typical 
"colorless object"? What is the lifetime of a typical "colorless object"? Do such objects 
have distinct electromagnetic structures? Are they hadron-like? Before more and better 
empirical facts about such objects become available, guesses and/or speculations may be 
helpful, provided that they agree with the existing data, and they are consistent with the 
fundamental theoretical knowledge — in particular consistent with the basic properties of the 
gluons (the direct gluon-gluon coupling prescribed by the QCD-Lagrangian, the confinement, 
and the non-conservation of gluon number, etc.). In this sense, we may wish to ask: Is it 
possible that the colorless objects are BTW-clusters which exist due to SOC in systems of 
interacting soft gluons? We are aware of the fact that the existence of SOC cannot (at least 
cannot yet) be derived from a basic theory such as QCD [Perhaps this can be and/or should 
be compared with the fact that the Gutenberg Richter law for earthquakes cannot (at least 
cannot yet) be derived from gravitational theory]. But, as in the case of earthquakes or 
any other open dynamical systems which leads to SOC, we can and we should ask: Can 
this be checked experimentally? Can this be done by looking for characteristic properties 
of SOC — in particular the SOC-fingerprints mentioned in Eqs.(|I|) and (0) in the relevant 
experiments? 

To answer these questions, it is useful to recall the following: Since the "colorless objects" 
are color- singlets which can exist inside and/or outside the proton, the interactions between 
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such color-singlets as well as those between such objects and "the mother proton" should 
be of Van der Waals type. Hence, it is expected that such a colorless object can be readily 
separated as an entire object from the mother proton in scattering processes in which the 
momentum-transfer is sufficient to overcome the binding energy due to the Van der Waals 
type of interactions. This means, in inelastic diffractive scattering the beam-particle (which 
is the virtual photon 7* in DIS) should have a chance to encounter one of the color-singlet 
gluon-clusters. For the reasons mentioned above, the struck colorless object can simply 
be "knocked out" and/or "carried away" by the beam-particle in such a collision event. 
Hence, it seems that the question whether "the colorless objects" are indeed BTW-clusters is 
something that can be answered experimentally. In this connection we recall that, according 
to the general theory of socSE, the size of a BTW-cluster is characterized by its dissipative 
energy, and in case of systems of interacting soft gluons associated with the proton, the 
dissipative energy carried by the BTW-cluster should be proportional to the energy fraction 
(xp) carried by the colorless object. Hence, if the colorless object can indeed be considered 
as a BTW-cluster due to SOC, we should be able to obtain information about the size- 
distribution of such color-singlet gluon-clusters by examing the xp-distributions of LRG 
events in the small-x^ region of DIS. 

Having this in mind, we now take a closer look at the measured! "diffractive structure 
function" F 2 D{3) (/3, Q 2 ; x P ) = J dtF? (A) (/3,Q 2 ; x P ,t). Here, we note that F 2 D{4) (/3, Q 2 ; x P , t) 
is relatedUM to the differential cross-section for large-rapidity-gap events 

dwdxpdt = -w ( " y + ^ 2 (A Q ; xp ' } ' (3) 

in analogy to the relationship between the corresponding quantities [namely d 2 a /(dxp dQ 2 ) 
and F2(xb, Q 2 )} for normal deep-inelastic electron-proton scattering events 

*' ^ 3 -» + £)*■.<«„,«•). (4) 



dx B dQ 2 x B Q 4 2 
The kinematical variables, in particular f3, Q 2 , x P and x p (in both cases) are directly 
measurable quantities, the definitions of which are shown in Fig.|I| together with the corre- 
sponding diagrams of the scattering processes. We note that, although these variables are 

8 



Lorentz-invariants, it is sometimes convenient to interpret them in a "fast moving frame", 
for example the electron-proton center-of-mass frame where the proton's 3-momentum P is 
large (i.e. its magnitude \P\ and thus the energy P° = (\P\ 2 + M 2 ) 1//2 is much larger than 
the proton mass M). While Q 2 characterizes the virtuality of the space-like photon 7*, xb 
can be interpreted, in such a "fast moving frame" (in the framework of the celebrated parton 
model), as the fraction of proton's energy P° (or longitudinal momentum \P\) carried by 
the struck charged constituent. 

We recall, in the framework of the parton model, ^(xg, Q 2 )/xb for "normal events" can 
be interpreted as the sum of the probability densities for the above-mentioned 7* to interact 
with a charged constituent of the proton. In analogy to this, the quantity F^ 3 \ft, Q 2 ; xp)/(3 
for LRG events can be interpreted as the sum of the probability densities for 7* to interact 
with a charged constituent which carries a fraction ft = xb/xp of the energy (or longitudinal 
momentum) of the colorless object, under the condition that the colorless object (which we 
associate with a system of interacting soft gluons) carries a fraction xp of proton's energy 
(or longitudinal momentum). We hereafter denote this charged-neutral and color- neutral 
gluon-system by Cq (in Regge pole models! this object is known as the "Pomeron"). Hence, 
by comparing Eq.(^) with Eq.(^) and by comparing the two diagrams shown in Fig.[TJa and 
Fig.|l]b, it is tempting to draw the following conclusions: 

The diffractive process is nothing else but a process in which the virtual photon 7* 
encounters a Cq, and ft is nothing else but the Bj or ken- variable with respect to Cg (this is 
why it is called xbc i n Ref.|[Llj|). This means, a diffractive e~p scattering event can be 
envisaged as an event in which the virtual photon 7* collides with "a Cg-target" instead of 
"the proton-target". Furthermore, since Cq is charge-neutral, and a photon can only directly 
interact with an object which has electric charges and/or magnetic moments, it is tempting 
to assign Cg an electromagnetic structure function F 2 c (/3, Q 2 ), and study the interactions 
between the virtual photon and the quark(s) and antiquark(s) inside Cg. In such a picture 
(which should be formally the same as that of Regge pole models!, if we would replace the 
Cg's by "pomerons") we are confronted with the following two questions: 
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First, is it possible and meaningful to discuss the xp-distributions of the Cq's without 
knowing the intrinsic properties, in particular the electromagnetic structures, of such ob- 
jects? 

Second, are the Cq's hadron-like, such that the electromagnetic structure of a (a typical, 
or an average) Cq can be studied in the same way as those for ordinary hadron? 

Since we wish to begin the quantitative discussion with something familiar to most of the 
readers in this community, and we wish to differentiate between the conventional-approach 
and the SOC-approach, we would like to discuss the second question here, and leave the 
first question to the next section. We recall (see in particular the last two papers in Ref.||), 
in order to see whether the second question can be answered in the affirmative, we need to 
know whether F^ 3 \[3, Q 2 ; x P ) can be factorized in the form 

F 2 D{3 \f3,Q 2 ;x P ) = f c (x P )F^,Q 2 ). (5) 

Here, f c {xp) plays the role of a "kinematical factor" associated with the "target Cg", and 
xp is the fraction of proton's energy (or longitudinal momentum) carried by Cq. [We could 
call f c (xp) "the Cq-Aux" — in exactly the same manner as in Regge pole models^, where it 
is called "the pomeron flux".] F£(f3,Q 2 ) is "the electromagnetic structure function of Cq" 
[the counterpart of F 2 (xb,Q 2 ) of the proton] which — in analogy to proton (or any other 
hadron) — can be expressed as 

= Y,e 2 [qmQ 2 )+mQ% (6) 

where fjf(<7f) stands for the probability density for 7* to interact with a quark (antiquark) 
of flavor i and electric charge e« which carries a fraction (3 of the energy (or longitudinal 
momentum) of Cq. It is clear that Eq.(^) should be valid for all xp- values in this kinematical 
region, that is, both the right- and the left-hand-side of Eq. (|6|) should be independent of the 
energy (momentum) carried by the "hadron" cj$. 

Hence, to find out experimentally whether the second question can be answered in the 
affirmative, we only need to check whether the data are in agreement with the assumption 
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that F£ (j3, Q 2 ) prescribed by Eqs.(|5D and (||) exists. For such a test, we take the existing 
data! and plot \og[F^ 3 \(3, Q 2 ; xp)/(3] against log/5 for different xp-values. We note, under 
the assumption that the factorization shown in Eq.(|^) is valid, the /3-dependence for a 
given Q 2 in such a plot should have exactly the same form as that in the corresponding 
logfF^/?, Q 2 )//3] vs log/5 plot; and that the latter is the analog of log[F 2 (xp, Q 2 )/xb] vs 
log xp plot for normal events. In Fig.||] we show the result of such plots for three fixed (Up- 
values (3.5, 20 and 65 GeV 2 , as representatives of three different ranges in Q 2 ). Our goal is 
to examine whether or how the /3-dependence of the function given in Eq.@ changes with 
xp. In principle, if there were enough data points, we should, and we could, do such a plot 
for the data-sets associated with every xp-value. But, unfortunately there are not so much 
data at present. What we can do, however, is to consider the /^-distributions in different 
xp-bins, and to vary the bin-size of xp, so that we can explicitly see whether/how the shapes 
of the /^-distributions change. The results are shown in FigPJ. The /^-distribution in the first 
row, corresponds to the integrated value i^P^, Q 2 ) shown in the literature^. Those in the 
second and in the third row are obtained by considering different bins and/or by varying 
the sizes of the bins. By joining the points associated with a given xp-interval in a plot 
for a given Q 2 , we obtain the /3-distribution for a Cq carrying approximately the amount 
of energy xpP°, encountered by a photon of virtuality Q 2 . Taken together with Eq.(f|) we 
can then extract the distributions q\ (/?, Q 2 ) and q c { {fi, Q 2 ) for this Q 2 -value, provided that 
-F 2 C (/?, Q 2 )//3 is independent of xp. But, as we can see in Fig.^|, the existing data! show 
that the xp-dependence of this function is far from being negligible! Note in particular 
that according to Eq.([|), by choosing a suitable fp{xp) we can shift the curves for different 
xp- values in the vertical direction (in this log-log plot); but we can never change the shapes 
of the (3 -distributions which are different for different xp- values! 

In order to see, and to realize, the meaning of the xp-dependence of the distributions 
of the charged constituents of Cq expressed in terms of -F 2 c (/5, Q 2 )/(3 in LRG events [see 
Eqs.(|J) and (^)], let us, for a moment, consider normal deep-inelastic scattering events 
in the xp-region where quarks dominate (xp > 0.1, say). Here we can plot the data for 
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log[i<2(xB, Q 2 )/xb] as a function of logxe obtained at different incident energies (P°'s) 
of the proton. Suppose we see, that at a given Q 2 , the data for xp-distributions taken 
at different values of P° are very much different. Would it still be possible to introduce 
F2(xb,Q 2 ) as "the electromagnetic structure function" of the proton, from which we can 
extract the x^-distribution of the quarks qi(xB,Q 2 ) at a given Q 2 1 The fact that it is not 
possible to assign an xp-independent structure function Q 2 )/ P to Cq which stands for 

the "Pomeron", and whose "flux" f c (xp) is expected to be independent of (3 and Q 2 , deserves 
to be taken seriously. It strongly suggests that the following picture cannot be true: "There 
exists an universal colorless object (call it Pomeron or or something else) the exchange 
of which describes diffractive scattering in general and DIS off proton in particular. This 
object is hadron-like in the sense that it has not only a typical size and a typical lifetime, 
but also a typical electromagnetic structure which can e.g. be measured and described by 
an electromagnetic structure function" . 

In summary of this section, we note that the empirical facts mentioned above show that 
no energy-independent electromagnetic structure function can be assigned to the expected 
universal colorless object Cq. This experimental fact is of considerable importance, because 
it is the first indication that, if there is an universal "colorless object", this object cannot be 
considered as an ordinary hadron. It has to be something else! In fact, as we shall see below, 
this property is closely related to the observation that such an object cannot have a typical 
size, or a typical lifetime. To be more precise, the fact that the data! cannot accommodate 
the simple factorization assumption shown in Eq.(|5|), in which a universal Pomeron flux 
with a unique hadron-like Pomeron structure function exists, can be considered as further 
support to the proposed SOC picture because a BTW-cluster, which has neither a typical size 
nor a typical lifetime, cannot have an universal static structure. With these characteristic 
properties of the colorless objects in mind, we see@ an overlap between the SOC-picture 
and the partonic picture for Pomeron and/or Pomeron and Regg eoni@ in which, beside the 
Pomeron, exchange of ( in general infinite number of) subleading trajectories are possible. In 
fact, it has been reportedi'S that very good agreement can be achieved between the datai! 
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and these types of models. Hence, in order to differentiate between the two approaches, 
it is not only useful but also necessary to examine the corresponding predictions for the 
dependence on the invariant momentum transfer t. This will be discussed in detail in 



Sections H - IT 



4. Distributions of the gluon-clusters 

After having seen that the existing data does not allow us to assign an energy-independent 
electromagnetic structure function to "the colorless object" such that the universal colorless 
object (eg) can be treated as an ordinary hadron, let us now come back to the first question 
in Section |3|, and try to find out whether it is never-the-less possible, and meaningful, to talk 
about the £p-distribution of Cq. As we shall see in this section, the answer to this question 
is Yes! Furthermore, we shall also see, in order to answer this question in the affirmative, 
we do not need the factorization mentioned in Eq.(|5|), and we do not need to know whether 
the Cq's are hadron-like. But, as we have already mentioned above, it is of considerable 
importance to discuss the second question so that we can understand the origin and the 
nature of the Cq's. 

In view of the fact that we do use the concept "distributions of gluons" in deep-inelastic 
lepton-hadron scattering, although the gluons do not directly interact with the virtual pho- 
tons, we shall try to introduce the notion "distribution of Cq" in a similar manner. In order 
to see what we should do for the introduction of such distributions, let us recall the following: 

For normal deep-inelastic e~p collision events, the structure function ^(^b, Q 2 ) can be 
expressed in terms of the distributions of partons, where the partons are not only quarks 
and antiquarks, but also gluons which can contribute to the structure function by quark- 
antiquark pair creation and annihilation. In fact, in order to satisfy energy-momentum- 
conservation (in the electron-proton system), the contribution of the gluons x g g(x g , Q 2 ) has 
to be taken into account in the energy-momentum sum rule for all measured Q 2 -values. Here, 
we denote by g(x g , Q 2 ) the probability density for the virtual photon 7* (with virtuality Q 2 ) 
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to meet a gluon which carries the energy (momentum) fraction x g of the proton, analogous 
to qi(xB,Q 2 ) [or qi(xB,Q 2 )] which stands for the probability density for this 7* to interact 
with a quark (or an antiquark) of flavor i and electric charge e, which carries the energy 
(momentum) fraction of the proton. We note, while both xb and x g stand for energy 
(or longitudinal momentum) fractions carried by partons, the former can be, but the latter 
cannot be directly measured. 

Having these, in particular the energy-momentum sum rule in mind, we immediately see 
the following: In a given kinematical region in which the contributions of only one category 
of partons (for example quarks for xb > 0.1 or gluons for xb < 10 -2 ) dominate, the structure 
function F 2 (xb, Q 2 ) can approximately be related to the distributions of that particular kind 
of partons in a very simply manner. In fact, the expressions below can be, and have been, 
interpreted as the probability-densities for the virtual photon 7* (with virtuality Q 2 ) to meet 
a quark or a gluon which carries the energy (momentum) fraction xb or x g respectively. 



The relationship between qi(xB,Q 2 ), g(x g ,Q 2 ) and F 2 (xb,Q 2 ) as they stand in Eq.(^) are 
general and formal (this is the case especially for that between g and F 2 ) in the following 
sense: Both qi(xs, Q 2 ) and g(x g , Q 2 ) contribute to the energy-momentum sum rule and both 
of them are in accordance with the assumption that partons of a given category (quarks or 
gluons) dominate a given kinematical region (here xb > 0.1 and xb < 10 -2 , respectively). 
But, neither the dynamics which leads to the observed Q 2 -dependence nor the relationship 
between x g and xb are given. This means, without further theoretical inputs, the simple 
expression for g(x g ,Q 2 ) as given by Eq.(^) is practically uselessl 

Having learned this, we now discuss what happens if we assume, in diffractive lepton- 
nucleon scattering, that the color-singlet gluon-clusters (cq's) dominate the small-x^ region 
(xb < 10 -2 , say). In this simple picture, we are assuming that the following is approximately 
true: The gluons in this region appear predominantly in form of Cq's. The interaction 
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or 



(7) 
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between the struck Cq and the rest of the proton can be neglected during the 7-Cq collision 
such that we can apply impulse-approximation to the Cq's in this kinematical region. That 
is, here we can introduce — in the same manner as we do for other partons (see Eq. ([?[)), 
a probability density Dg(xp\/3, Q 2 ) for 7* in the diffractive scattering process to "meet" a 
Cq which carries the fraction xp of the proton's energy P° = (\P\ 2 + M 2 ) 1 / 2 « \P\ (where 
P is the momentum and M is the mass of the proton). In other words, in diffractive 
scattering events for processes in the kinematical region xp < 10~ 2 , we should have, instead 
of g(x g ,Q 2 ), the following: 

F^\^Q 2 -x P ) 



x P 

?0 



D s (xp\(3,Q 2 ). (8) 



Here, xpP is the energy carried by Cg, and j3 indicates the corresponding fraction carried by 
the struck charged constituent in Cg. In connection with the similarities and the differences 
between qi(xp,Q 2 ), g(xp,Q 2 ) in Eq.(^) and D s (xp\[3,Q 2 ) in Eq.(|J), it is useful to note 
in particular the significant difference between x g and xp, and thus that between the x g - 
distribution g(x g ,Q 2 ) of the gluons and the xp-distribution Ds(xp\(3, Q 2 ) of the Cq's: Both 
x g and xp are energy (longitudinal momentum) fractions of charge- neutral objects, with 
which 7* cannot directly interact. But, in contrast to x g , xp can be directly measured in 
experiments, namely by making use of the kinematical relation 

Q 2 + M * /n, 

xp K WTWf 0) 

and by measuring the quantities Q 2 , M 2 and W 2 in every collision event. Here, Q, M x and 
W stand respectively for the invariant momentum-transfer from the incident electron, the 
invariant-mass of the final hadronic state after the 7* — Cq collision, and the invariant mass 
of the entire hadronic system in the collision between 7* and the proton. Note that xp = 
/3xp, hence f3 is also measurable. This means, in sharp contrast to g(x g ,Q 2 ), experimental 
information on Ds(xp\/3, Q 2 ) in particular its xp-dependence can be obtained — without 
further theoretical inputs] 
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5. The first SOC-fingerprint: Spatial scaling 

We mentioned at the beginning of Section |3], that in order to find out whether the concept 
of SOC indeed plays a role in diffractive DIS we need to check the fingerprints of SOC 
shown in Section |2|, and that such tests can be made by examing the corresponding cluster- 
distributions obtained from experimental data. We are now ready to do this, because we 
have learned in Sections |3] and |], that it is not only meaningful but also possible to extract 
xp-distributions from the measured diffractive structure functions, although the Cq's cannot 
be treated as hadrons. In fact, as we can explicitly see in Eqs.(]^) and (Q), in order to 
extract the xp-dependence of the Cq's from the data, detailed knowledge about the intrinsic 
structure of the Cq's are not necessary. 

Having these in mind, we now consider as a function of xp for given values of 
(3 and Q 2 , and plot F^ 3 \(3,Q 2 ;xp)/xp against xp for different sets of (3 and Q 2 . The 
results of such log-log plots are shown in Fig|| As we can see, the dataifl suggest that the 
probability-density for the virtual photon 7* to meet a color-neutral and charged-neutral 
object Cq with energy (longitudinal momentum) fraction xp has a power-law behavior in 
xp, and the exponent of this power- law depends very little on Q 2 and (3. This is to be 
compared with D S (S) in Eq.( [[]), where S, the dissipative energy (the size of the BTW- 
cluster) corresponds to the energy of the system Cq. The latter is xpP°, where P° is the 
total energy of the proton. 

It means, the existing data!'! show that Ds(xp\f3, Q 2 ) exhibits the same kind of power- 
law behavior as the size-distribution of BTW-clusters. This result is in accordance with the 
expectation that self-organized critical phenomena may exist in systems of interacting soft 
gluons in diffractive deep-inelastic electron-proton scattering processes. 

We note, up to now, we have only argued (in Section P that such gluon-systems are open, 
dynamical, complex systems in which SOC may occur, and we have mentioned (in Section 0) 
the ubiquitousness of SOC in Nature. Having seen the experimental evidence that one of the 
"SOC-fingerprints" (which are necessary conditions for the existence of SOC) indeed exists, 
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let us now take a second look at such gluon-systems from a theoretical standpoint. Viewed 
from a "fast moving frame" which can for example be the electron-proton c.m.s. frame, such 
systems of interacting soft gluons are part of the proton (although color- singlets can also be 
outside the confinement region). Soft gluons can be intermittently emitted or absorbed by 
gluons in such a system, as well as by gluons, quarks and antiquarks outside the system. The 
emission- and absorption-processes are due to local interactions prescribed by the well-known 
QCD-Lagrangian (here "the running coupling constants" are in general large, because the 
distances between the interacting colored objects cannot be considered as "short"; remember 
that the spatial dimension of be much larger than that of a hadron!). In this 

connection, it is useful to keep the following in mind: Due to the complexity of the system, 
details about the local interactions may be relatively unimportant, while general and/or 
global features — for example energy-flow between different parts (neighbors and neighbor's 
neighbors . . .) of the system — may play an important role. 

How far can one go in neglecting dynamical details when one deals with such open com- 
plex systems? In order to see this, let us recall how Bak and Sneppen0 succeeded in modeling 
some of the essential aspects of The Evolution in Nature. They consider the "fitness" of 
different "species", related to one another through a "food chain", and assumed that the 
species with the lowest fitness is most likely to disappear or mutate at the next time-step 
in their computer simulations. The crucial step in their simulations that drives evolution 
is the adaption of the individual species to its present environment (neighborhood) through 
mutation and selection of a fitter variant. Other interacting species form part of the envi- 
ronment. This means, the neighbors will be influenced by every time-step. The result these 
authors obtained strongly suggests that the process of evolution is a self-organized critical 
phenomenon. One of the essential simplifications they made in their evolution mo 

delsSH 

is 

the following: Instead of the explicit connection between the fitness and the configuration of 
the genetic codes, they use random numbers for the fitness of the species. Furthermore, as 
they have pointed out in their papers, they could in principle have chosen to model evolution 
on a less coarse-grained scale by considering mutations at the individual level rather than 
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on the level of species, but that would make the computation prohibitively difficult. 

Having these in mind, we are naturally led to the questions: Can we consider the cre- 
ation and annihilation processes of colorless systems of interacting soft gluons associated 
with a proton as "evolution" in a microscopic world? Before we try to build models for a 
quantitative description of the data, can we simply apply the existing evolution modelsSll 
to such open, dynamical, complex systems of interacting soft gluons, and check whether 
some of the essential features of such systems can be reproduced? 

To answer these questions, we now report on the result of our first trial in this direction: 
Based on the fact that we know very little about the detailed reaction mechanisms in such 
gluon-systems and practically nothing about their structures, we simply ignore them, and 
assume that they are self-similar in space (this means, color-singlet gluon-clusters (cq) can 
be considered as clusters of Cq's and so on). Next, we divide them into an arbitrary given 
number of subsystems Sj (which may or may not have the same size). Such a system is open, 
in the sense that neither its energy nor its gluon-number n t has a fixed value. Since we do 
not know, in particular, how large the £j's are, we use random numbers. As far as the n(s 
are concerned, since we do not know how these numbers are associated with the energies in 
the subsystems Sj, except that they are not conserved quantities, we just ignore them, and 
consider only the Si's. As in Ref.[|T7| or in Ref.|l8|, the random number of this subsystem 



as well as those of the fixedO or random (see the first paper of Ref.|TjJ) neighbors will be 



changed at every time-step. Note, this is how we simulate the processes of energy flow due to 
exchange of gluons between the subsystems, as well as those with gluons/quarks/antiquarks 
outside the system. In other words, in the spirit of Bak and SneppenEll we are neglecting 
the dynamical details totally. Having in mind that, in such systems, the gluons as well as 
the subsystems (sj's) of gluons are virtual (space-like), we can ask: "How long can such a 
colorless subsystem s« of interacting soft gluons exist, which carries energy 8(1" According to 
the uncertainty principle, the answer should be: "The time interval in which the subsystem 
Si can exist is proportional to and this quantity can be considered as the lifetime 

Tj of Si" In this sense, those colorless subsystems of gluons are expected to have larger 



probabilities to mutate when they are associated with higher energies, and thus shorter 
"lifetimes". Note that the basic local interaction in this self-organized evolution process 
is the emission (or absorption) of gluons by gluons prescribed by the QCD-Lagrangian - 
although the detailed mechanisms (which can in principle be explicitly written down by 
using the QCD-Lagrangian) do not play a significant role. 

In terms of the evolution model0lll we now call Sj the "species" and identify the corre- 
sponding lifetime 73 as the "fitness of s" . Because of the one-to-one correspondence between 
Ti and Ei, where the latter is a random number, we can also directly assign random numbers 
to the Tj's instead. From now we can adopt the evolution mo and note that, at the 

start of such a process (a simulation), the fitness on average grow, because the least fit are 
always eliminated. Eventually the fitness do not grow any further on average. All gluons 
have a fitness above some threshold. At the next step, the least fit species (i.e. the most 
energetic subsystem Sj of interacting soft gluons), which would be right at the threshold, 
will be "replaced" and starts an avalanche (or punctuation of mutation events), which is 
causally connected with this triggering "replacement". After a while, the avalanche will 
stop, when all the fitnesses again will be over that threshold. In this sense, the evolution 
goes on, and on, and on. As in Refs.[|17| and [IT3], we can monitor the duration of every 



avalanche, that is the total number of mutation events in everyone of them, and count how 
many avalanches of each size are observed. The avalanches mentioned here are special cases 
of those discussed in Section |[ Their size- and lifetime-distributions are given by Eq.(P 
and Eq.(p]), respectively. Note in particular that the avalanches in the Bak-Sneppen model 
correspond to sets of subsystems Sj, the energies (e*) of which are too high "to be fit for the 
colorless systems of low-energy gluons". It means, in the proposed picture, what the virtual 
photon in deep-inelastic electron-proton scattering "meet" are those "less fit" one — those 
who carry "too much" energy. In a geometrical picture this means, it is more probable for 
such "relatively energetic" color-singlet gluon-clusters (cq) to be spatially further away from 
the (confinement region of) the proton. 

There exists, in the mean time, already several versions of evolution mo delsS@ based 



on 
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the original idea of Bak and SneppenO. Although SOC phenomena have been observed in all 
these cases the slopes of the power-law distributions for the avalanches are different 

in different models — depending on the rules applied to the mutations. The values range 
from approximately —1 to approximately —2. Furthermore, these models000 seem to 
show that neither the size nor the dimension of the system used for the computer simulation 
plays a significant role. 

Hence, if we identify the colorless charge-neutral object encountered by the virtual 
photon 7* with such an avalanche, we are identifying the lifetime of Cg with T, and the "size" 
(that is the total amount of dissipative energy in this "avalanche") with the total amount 
of energy of cj$. Note that the latter is nothing else but xpP°, where P° is the total energy 
of the proton. This is how and why the ^-distribution in Eq. ([]]) and the xp-distribution of 
Ds(%p\/3, Q 2 ) in Eq.(||) are related to each other. 



6. The second fingerprint: Temporal scaling 

In this section we discuss in more detail the effects associated with the time-degree-of- 
freedom. In this connection, some of the concepts and methods related to the two questions 
raised in Section ||| are of great interest. In particular, one may wish to know why the parton- 
picture does not work equally well for hadrons and for Cq's. The answer is very simple: The 
time-degree of freedom cannot be ignored when we apply impulse-approximation, and the 
applicability of the latter is the basis of the parton-model. We recall that, when we apply 
the parton model to stable hadrons, the quarks, antiquarks and gluons are considered as free 
and stable objects, while the virtual photon 7* is associated with a given interaction-time 
T~int{Q 2 , xb) characterized by the values Q 2 and xb of such scattering processes. We note 
however that, this is possible only when the interaction-time r- int is much shorter than the 
corresponding time-scales related to hadron-structure (in particular the average propagation- 
time of color-interactions in hadrons). Having these in mind, we see that, we are confronted 
with the following questions when we deal with Cq's which have finite lifetimes: Can we 
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consider the Cq's as "/ree" and "stable" particles when their lifetimes are shorter than the 
interaction-time r int (Q 2 , Can we say that a 7*-Cq collision process takes place, in which 
the incident 7* is absorbed by one or a system of the charged constituents of Cq, when the 
lifetime T of Cq is shorter than r int (Q 2 , £b)? 

Since the notion "stable objects" or "unstable objects" depends on the scale which is 
used in the measurement, the question whether be considered parton (in the 

sense that it can be considered as a free "stable object" during the 7*-Cq interaction) depends 
very much on the interaction-time Ti nt (Q 2 ', x b) ■ Here, for given values of Q 2 , xb, and thus 
T int(Q 2 ,XB), only those Cq's whose lifetimes (T's) are greater than T int (Q 2 ,XB) can absorb 
the corresponding 7*. That is to say, when we consider diffractive electron-proton scattering 
in kinematical regions in which Cq's dominate, we must keep in mind that the measured 
cross-sections (and thus the diffractive structure function i 7 ^ 3 -*) only include contributions 
from collision-events in which the condition T > Ti n t(Q 2 ,XB) is satisfied! 

We note that r int can be estimated by making use of the uncertainty principle. In fact, 
by calculating l/q° in the above-mentioned reference frame, we obtain 

A\P\xb , 1n * 

Tint = q^T^? (10) 

which implies that, for given \P\ and Q 2 values, 

ri nt oc x B , for x B < 1. (11) 

This means, for diffractive e~p scattering events in the small-x^ region at given \P\ and 
Q 2 values, xb is directly proportional to the interaction time T mt . Taken together with the 
relationship between r int and the minimum lifetime T(min) of the Cq's mentioned above, 
we reach the following conclusion: The distribution of this minimum value, T(min) of the 
Cq's which dominate the small-x^ (xb < 1CT 2 , say) region can be obtained by examining 
the xe-dependence of F^ 3 \f3,Q 2 ; Xp)/ (3 discussed in Eqs.(^), (^) and in Fig.^. This is 
because, due to the fact that this function is proportional to the quark (antiquark) distribu- 
tions q1{qi c ) which can be directly probed by the incident virtual photon 7*, by measuring 
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F^ 3 \f3, Q 2 , xp)/ (3 as a function of xb = (3xp, we are in fact asking the following question: 
Do the distributions of the charged constituents of Cq depend on the interaction time r int , 
and thus on the minimum lifetime T(min) of the to be detected Cq? We use the identity 
xb = (3xp and plot the quantity {j3, Q 2 ; xp)/(3 against the variable xb for fixed values 
of (3 and Q 2 . The result of such a log-log plot is given in Fig.^. It shows not only how 
the dependence on the time-degree-of- freedom can be extracted from the existing datai, but 
also that, for all the measured values of (3 and Q 2 , the quantity 

P (M), F ^f^ (12, 

is approximately independent of (3, and independent of Q 2 . This strongly suggests that the 
quantity given in Eq.([E|) is associated with some global features of Cq — consistent with the 
observation made in Section |] which shows that it cannot be used to describe the structure 
of Cq. This piece of empirical fact can be expressed by setting p(xb\(3, Q 2 ) ~ p(xb)- By 
taking a closer look at this log-log plot, as well as the corresponding plots for different sets 
of fixed (3- and Q 2 - values (such plots are not shown here, they are similar to those in Fig.^J), 
we see that they are straight lines indicating that p(xs) obeys a power-law. What does this 
piece of experimental fact tell us? What can we learn from the distribution of the lower 
limit of the lifetimes (of the gluon-systems Cg's)? 

In order to answer these questions, let us, for a moment, assume that we know the 
lifetime-distribution Dt(T) of the Cq's. In such a case, we can readily evaluate the integral 



J[r int (x B )] = / D T (T)dT, (13) 

and thus obtain the number density of all those clusters which live longer than the interaction 
time r int (xB). Hence, under the statistical assumption that the chance for a 7* to be absorbed 
by one of those Cq's of lifetime T is proportional to Dt(T) (provided that T\ nt {Q 2 , x b) < T, 
otherwise this chance is zero), we can then interpret the integral in Eq.fJTS]) as follows: 
[fint(Q 2 , xb)} oc p(xb) is the probability density for 7* [associated with the interaction-time 
Tint{%B)} to be absorbed by Cq's. Hence, 
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D T (x B ) oc -^~p( x b)- (14) 

This means in particular, the fact that p(xs) [introduced in Eq. (|T2|) l obeys a power- law in 
Xb implies that Dx{T) obeys a power-law in T. Such a behavior is similar to that shown 
in Eq.( §). In order to see the quality of this power-law behavior of Dt, and the quality of 
its independence of Q 2 and (3, we compare the above-mentioned behavior with the existing 
dataiH. In Fig.^, we show the log-log plots of d/dxB[p(xB)} against xb- In doing this plot, 
we keep the definition of p(xs\P, Q 2 ) given in Eq.(|T2|) and its weak (3- and Q 2 -dependence in 
mind, and we note that d/dxB[p(xs)] is approximately F^ 3 \[3 } Q 2 ;xb/(3)/((3xb), provided 
that p(xs\f3, Q 2 ) shows a power-law behavior in xb- Here, we not only see that the quality 
of the power-law behavior of in T is intimately related to the quality of the power- 
law behavior of F^ 3 \(3,Q 2 ;xb/(3)/((3xb) in but also how weak the (3- and the Q 2 - 
dependence are. 



7. SOC-fingerprints in inelastic diffractive "y*p—,jp—,pp—, and pp— scattering 
processes 

We have seen, in Sections |5] and |6|, that in diffractive deep-inelastic electron-proton scat- 
tering, the size- and the lifetime-distributions of the color-singlet gluon-clusters (eg) obey 
power- laws, and that the exponents depend very little on the variables j3 and Q 2 . We inter- 
preted the power-law behaviors as the fingerprints of SOC which are expected to manifest 
themselves in systems of interacting soft gluons (which play the dominating role in diffrac- 
tive DIS). This expectation is based on the fact that the fundamental properties of gluons 
(in particular the direct gluon-gluon coupling, the confinement, and the non-conservation of 
gluon-numbers) show that the necessary conditions for the existence of SOC in systems of 
interacting gluons are fulfilled, and the fact that (as we can see in various open dynamical 
complex systems) power-law behavior in size- and lifetime-distributions are indeed reliable 
indicators for the existence of SOC. In this sense, the existence of such power-law behavior 
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can be understood in terms of the QCD-based SOC-picture, although it is not (at least not 
yet) possible to derive the power-law behavior of the size- and lifetime-distributions, and to 
calculate the exponents by using non-perturbative QCD. But can the observed approximate 
independence (or weak dependence) of the exponents on Q 2 and (3 also be understood in 
terms of the QCD-based SOC-picture mentioned above? In particular, what do we expect to 
see in photoproduction processes where the associated value for Q 2 is approximately zero? 
We note that the possible relationship between the Q 2 > a few GeV 2 case and the Q 2 ~ 
case in diffractive scattering is of considerable interest for many reasons. One of them is the 
fact that, by comparing these two cases, we can see the fundamental difference between the 
conventional (pQCD-corrected parton model plus Regge phenomenology) picture and the 
proposed QCD-based SOC-picture for inelastic diffractive scattering. In the conventional 
picture, the Q 2 > a few GeV 2 case is "hard" and thus should be described by concepts and 
methods of parton model and pQCD, while the Q 2 ~ case is "soft" and thus should be 
understood in terms of Regge poles. What are the predictions of the proposed SOC-picture? 
What do the experimental data tell us in this connection? Would a systematic comparison 
of the existing data at different Q 2 -values — including those near Q 2 = be useful in un- 
derstanding the underlying reaction mechanism(s) of diffractive scattering in general, and 
differentiate between the conventional and the proposed picture in particular? 

In order to answer these questions, let us recall the space-time aspects of the collision 
processes which are closely related to the above-mentioned power-law behaviors. Viewed in 
a fast moving frame (e.g. the c.m.s. of the colliding electron and proton), the states of the 
interacting soft gluons originating from the proton are self-organized. The Cq caused by local 
perturbations and developed through "domino effects" are BTW-avalanches (see Sections [I] 
and |Sp, the size-distribution of which [see Eqs.(^) and ([]])] are given by Fig.|3|. This explicitly 
shows that there are Cq's of all sizes, because a power-law size-distribution implies that there 
is no scale in size. Recall that, since such Cq's are color-singlets, their spatial extensions can 
be much larger than that of the proton, and thus they can be "seen" also outside the proton 
by a virtual photon originating from the electron. In other words, what the virtual photon 
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encounters is a cloud of Cq's, everyone of which is in general partly inside and partly outside 
the proton. 

The virtual photon, when it encounters a Cq, will be absorbed by the charged constituents 
(quarks and antiquarks due to fluctuations of the gluons) of the gluon-system. Here it 
is useful to recall that in such a space-time picture, Q 2 is inversely proportional to the 
transverse size, and xb is a measure of the interaction time [See Eqs. (|I]) and (|TT1) in 
Section |(| of the virtual photon. It is conceivable, that the values for the cross-sections for 
virtual photons (associated with a given Q 2 and a given xb) to collide with Cq's (of a given 
size and a given lifetime) may depend on these variables. But, since the processes of self- 
organization (which produce such Cq's) take place independent of the virtual photon (which 
originates from the incident electron and enters "the cloud" to look for suitable partners), 
the power-law behaviors of the size- and lifetime-distributions of the Cq's are expected to 
be independent of the properties associated with the virtual photon. This means, by using 
7*'s associated with different values of Q 2 to detect Cq's of various sizes, we are moving up 
or down on the straight lines in the log-log plots for the size- and lifetime distributions, the 
slopes of which do not change. In other words, the observed approximative Q 2 -independence 
of the slope in the above-mentioned log-log plots of the data can be considered as a natural 
consequence of the QCD-based SOC picture. 

As far as the /^-dependence is concerned, we recall the results obtained in Sections |3| 
and (|, which explicitly show the following: The Cq's cannot be considered as hadrons. In 
particular, it is neither possible nor meaningful to talk about "the electromagnetic structure 
of the (or a typical, or an average) Cq". This is not only because the power-law behavior 
of the size and the lifetime distributions of such Cq's implies that such objects — although 
they are color- singlets — can have neither a typical (an average) size, nor a typical (an 
average) lifetime, but also because of the following fact: When F^ 3 \f3, Q 2 ; Xp)/j3 which is 
usually knowni as "the electromagnetic structure function of the colorless object exchanged 
in diffractive scattering processes" is plotted as functions of (3 (cf. Fig.|^), we see a rather 
significant xp-dependence. This does not mean, however, that the measured /3-dependence 
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of F®^ cannot provide us with any further information on the electromagnetic properties 
of the color-singlet gluon-clusters (cj)! This is because, the Cq's which play the dominating 
role in diffractive scattering are color-singlets, hence, even when such clusters are BTW- 
avalanches which have neither a typical size nor a typical lifetime, a set of such clusters 
with given size- and lifetime-distributions can nevertheless be considered as a specific set 
of hadrons with distinct properties. Hence, when the /3-dependence of F^ 3 \(3, Q 2 ; xp)/(3 is 
examined in inelastic diffractive scattering processes, the electromagnetic properties of such 
a set of hadrons are probed by the incident (virtual or real, depending on the Q 2 -value of 
the event) photons. In this connection, it is perhaps useful to consider the /^-distribution 
integrated over xp 

For the purpose of comparing SOC-fmgerprints obtained at different Q 2 -values, we are 
interested much more in measurable quantities in which the integrations over (3 have been 

carried out. A suitable candidate for this purpose is the differential cross-section 

1 d 2 a D 
xp dQ 2 dxp 



(3Q 4 \ 2 J x P 

dP J¥ Ds{Xp1 ^ Q2) (15) 

Together with Eqs.(|3]) and (^), we see that this cross-section is nothing else but the effective 
/3-weighted xp-distribution Ds(xp\Q 2 , (3) of the Cq's. Note that the weighting factors shown 
on the right-hand-side of Eq.(|l5|) are simply results of QED! Next, we use the dat all for 
F®^ which are available at present, to do a log- log plot for the integrand of the expression 



in Eq.([l5|) as a function of xp for different values of (3 and Q 2 . This is shown in in Fig.^a 
and Fig.||c. Since the absolute values of this quantity depend very much, but the slope of 
the curves very little on (3, we carry out the integration as follows: We first fit every set 
of the data separately. Having obtained the slopes and the intersection points, we use the 
obtained fits to perform the integration over (3. The results are shown in the 

. / 1 d 2 a D \ 
log — — versus log (xp) 



xp dQ 2 dxp 
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plots of Fig.|6|b. These results show the Q 2 -dependence of the slopes is practically negligible, 
and that the slope is approximately —1.95 for all values of Q 2 . 

Furthermore, in order to see whether the quantity introduced in Eq.(|l5D is indeed useful, 
and in order to perform a decisive test of the Q 2 -independence of the slope in the power-law 
behavior of the above-mentioned size-distributions, we now compare the results in deep- 
inelastic scatteringi'S with those obtained in photoproduction0, where LRG events have also 
been observed. This means, as in diffractive deep-inelastic scattering, we again associate 
the observed effects with colorless objects which are interpreted as a system of interacting 
soft gluons originating from the proton. In order to find out whether it is the same kind of 
gluon-clusters as in deep-inelastic scattering, and whether they "look" very much different 
when we probe them with real (Q 2 = 0) photons, we replot the existing da/dM 2 data@ for 
photoproduction experiments performed at different total energies, and note the kinematical 

relationship between M 2 , W 2 and x P (for Q 2 < M 2 and \t\ < Mj): 
M 2 

**»w (16) 

The result of the corresponding 

log ( j—) versus log (x P ) 

\xp dxp J 

plot is shown in Fig.^. The slope obtained from a least-square fit to the existing data@ is 
-1.98 ±0.07. 

The results obtained in diffractive deep-inelastic electron-proton scattering and that for 
diffractive photoproduction strongly suggest the following: The formation processes of Cq's 
in the proton is due to self-organized criticality, and thus the spatial distributions of such 
clusters — represented by the xp-distribution — obey power-laws. The exponents of such 
power-laws are independent of Q 2 . Since 1/Q 2 can be interpreted as a measure for the 
transverse size of the incident virtual photon, the observed Q 2 -independence of the exponents 
can be considered as further evidence for SOC — in the sense that the self-organized gluon- 
cluster formation processes take place independent of the photon which is "sent in" to detect 
the clusters. 

27 



Having these results, and the close relationship between real photons and hadrons in 
mind, we are immediately led to the following questions: What shall we see, when we 
replace the (virtual or real) photon by a hadron — a proton or an antiproton? (See in 
this connection Fig.|], for the notations and the kinematical relations for the description of 
such scattering processes.) Should we not see similar behaviors, if SOC in gluon-systems is 
indeed the reason for the occurrence of Cq's which can be probed experimentally in inelastic 
diffractive scattering processes? To answer these questions, we took a closer look at the 
available single diffractive proton-proton and proton-antiproton scattering dataiil; and in 
order to make quantitative comparisons, we plot the quantities which correspond to those 
shown in Fig.||b and Fig.[7[ These plots are shown in Fig.|9|a and Fig.pD. In Fig.|9|a, we see 
the double differential cross-section (l/xp)d 2 a/(dtdxp) at four different t-values. In Fig.||b, 
we see the integrated differential cross-section (1/ 'xp)da / 'dxp. Note that, here 

x P « M 2 Js, (17) 

where a/s is the total c.m.s. energy of the colliding proton-proton or antiproton-proton 
system. Here, the integrations of the double differential cross-section over t are in the 
ranges in which the corresponding experiments have been performed. (The extremely weak 
energy-dependence has been ignored in the integration.) The dashed lines in all the plots 
in Figs.|9|a and |9|b stand for the slope —1.97 which is the average of the slope obtained from 
the plots shown in Figs.^o and |7|. This means, the result shows exactly what we expect 
to see: The fingerprints of SOC can be clearly seen also in proton- and antiproton-induced 
inelastic diffractive scattering processes, showing that such characteristic features are indeed 
universal and robust! 

We are thus led to the following conclusions. Color-singlet gluon-clusters (cq) can be 
formed in hadrons as a consequence of self-organized criticality (SOC) in systems of in- 
teracting soft gluons. In other words, "the colorless objects" which dominate the inelastic 
diffractive scattering processes are BTW-avalanches (BTW-clusters). Such Cq's are in general 
distributed partly inside and partly outside the confinement region of the "mother-hadron" . 
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Since the interactions between the Cg's and other color-singlet objects (including the target 
proton) should be of Van der Waals type, it is expected that such an object can be readily 
driven out of the above-mentioned confinement region by the beam-particle in geometrically 
more peripheral collisions. This has been checked by examing inelastic single-diffractive 
scattering processes at high energies in which virtual photon, real photon, proton, and an- 
tiproton are used as beam particles. The result of this systematic check shows that the 
universal distributions of such Cq's can be directly extracted from the data. In particular, 
the fact that xp is the energy fraction carried by the struck Cq's, and the fact that the 
xp-distributions are universal, it is tempting to regard such xp-distributions as "the parton- 
distributions" for diffractive scattering processes. Having seen this, it is also tempting to 
ask: Can we make use of such "parton-distributions" to describe and/or to predict the mea- 
surable cross-sections in inelastic diffractive scattering processes? This and other related 
questions will be discussed in the following sections. 

8. Diffractive scattering in high-energy collisions and diffraction in optics 

It might sound strange, but it is true that people working in this field of physics often do 
not agree with one another on the question: What is "Diffractive Scattering in High- Energy 
Collisions"? In this paper, we have, until now, simply adopted the currently! popular defini- 
tion of "inelastic diffractive scattering processes" . That is, when we talked about "inelastic 
diffractive scattering" we were always referring to processes in which "colorless objects" are 
"exchanged" . In other words, until now, the following question has not been asked: Are the 
above-mentioned "inelastic diffractive scattering processes" indeed comparable with diffrac- 
tion in Optics, in the sense that the beam particles should be considered as waves, and 
the target-proton together with the associated (in whatever manner) colorless objects can 
indeed be viewed "scattering screen"? 

This question will be discussed in the present and the subsequent sections, together with 
the existing; data!! for the double differential cross-section d 2 a/dt d(M%/ s) for proton- 
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proton and antiproton-proton collisions (where t is the 4-momentum-transfer squared, M x is 
the missing-mass, and yfs is the total c.m.s. energy). The purpose of this investigation is to 
find out: "Can the observed t-dependence and the (M 2 /s)-dependence of d? 'a / ' dt d(M 2 / s) in 
the given kinematic range (0.2 GeV 2 <|i|< 3.25 GeV 2 , lQGeV<y/s< 630 GeV, and M 2 / ' s< 
0.1) be understood in terms of the well-known concept of diffraction in Optics? The answer 
to this question is of particular interest for several reasons: 

(a) High-energy proton-proton and proton-antiproton scattering at small momentum 
transfer has played, and is still playing a very special role in understanding diffraction 
and/or diffractive dissociation in lepton-, photon- and hadr on- induced reaction 

Many experiments have been performed at various incident energies for elastic and inelas- 
tic diffractive scattering processes. It is known that the double differential cross section 
d 2 a / dt d(M 2 / s) is a quantity which can yield much information on the reaction mecha- 
nism^) and/or on the structure of the participating colliding objects. In the past, the 
t-, M x - and s-dependence of the differential cross-sections for inelastic diffractive scatter- 
ing processes has been presented in different forms, where a number of interesting features 
have been observedfHS For example, it is seen that, the t-dependence of d 2 o /dtdM 2 at 
fixed s depends very much on My. the M 2 -dependence of d 2 a /dtdM 2 at fixed t depends 
on s. But, when d 2 a j dt d[M 2 j s) is plotted as function of M 2 /s at given t -values (in 
the range 0.2 GeV 2 < \t\ < 3.25 GeV 2 ) they are approximately independent of s\ What 
do these observed striking features tell us? The first precision measurement of this quan- 
tity was published more than twenty years ago!. Can this, as well as the more recent 
d 2 a/dt d(M 2 / s)-d&taM be understood theoretically? 

(b) The idea of using optical and/or geometrical analogies to describe high-energy 
hadron-nucleus and hadron-hadron collisions at small scattering angles has been discussed 
by many authors@S many years ago. It is shown in particular that this approach is very 
useful and successful in describing elastic scattering. However, it seems that, until now, no 
attempt has been made to describe the datai@ by performing quantitative calculations for 
d 2 a/dtd(M 2 /s) by using optical geometrical analogies. It seems worthwhile to make such an 
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attempt. This is because, it has been pointed outO very recently, that the above-mentioned 
analogy can be made to understand the observed ^-dependence in da/dt. 

(c) Inelastic diffractive pp- and pp-scattering belongs to those soft processes which have 
also been extensivly discussed in the well-known Regge-pole approachi'iil. The basic idea 
of this approach is that colorless objects in form of Regge trajectories (Pomerons, Reggeons 
etc.) are exchanged during the collision, and such trajectories are responsible for the dy- 
namics of the scattering processes. In this approach, it is the i-dependence of the Regge 
trajectories, the t-dependence of the corresponding Regge residue functions, the properties 
of the coupling of the contributing trajectories (e.g. triple Pomeron or Pomeron-Reggeon- 
Pomeron coupling), and the number of contributing Regge trajectories which determine the 
experimentally observed t- and M^-dependence of d 2 a/dtd(M 2 /s). A number of Regge-pole 
models§! have been proposed, and there exist good fitJi to the data. What remains to 
be understood in this approach is the dynamical origin of the Regge-trajectories on the one 
hand, and the physical meaning of the unknown functions (for example the t-dependence of 
any one of the Regge-residue functions) on the other. It has been pointed outBS, that there 
may be an overlap between the "Partons in Pomeron and Reggeons" picture and the SOC- 
picture0, and that one way to study the possible relationship between the two approaches 
is to take a closer look at the double differential cross-section d 2 a/dtd(M 2 /s). 



9. Optical diffraction off dynamical complex systems 

Let us begin our discussion on the above-mentioned questions by recalling that the concept 
of "diffraction" or "diffractive scattering" has its origin in Optics, and Optics is part of 
Electrodynamics, which is not only the classical limit, but also the basis of Quantum Elec- 
trodynamics (QED). Here, it is useful to recall in particular the following: Optical diffraction 
is associated with departure from geometrical optics caused by the finite wavelength of light. 
Frauenhofer diffraction can be observed by placing a scatterer (which can in general be a 
scattering screen with more than one aperture or a system of scattering objects) in the 



31 



path of propagation of light (the wavelength of which is less than the linear dimension of 
the scatterer) where not only the light-source, but also the detecting device, are very far 
away from the scatterer. The parallel incident light-rays can be considered as plane waves 
(characterized by a set of constants k, w = \k\, and u say, which denote the wave vector, 
the frequency and the amplitude of a component of the electromagnetic field respectively in 
the laboratory frame). After the scattering, the scattered field can be written in accordance 
with Huygens' principle as 

p i\k'\R 

u P = —f{k,k'). (18) 

Here, up stands for a component of the field originating from the scatterer, k! is the wave 
vector of the scattered light in the direction of observation, \k'\ = u' is the corresponding 
frequency, R is the distance between the scatterer and the observation point P, and f(k, k') 
is the (unnormalized) scattering amplitude which describes the change of the wave vector in 
the scattering process. By choosing a coordinate system in which the z-axis coincides with 
the incident wave vector k, the scattering amplitude can be expressed as followJ^Eii 

f(q) = j^y 2 JJd 2 ba(b)e-^ b . (19) 

Here, q = k' — k determines the change in wave vector due to diffraction; b is the impact 
parameter which indicates the position of an infinitesimal surface element on the wave-front 
"immediately behind the scatterer" where the incident wave would reach in the limit of 
geometrical optics, and a{b) is the corresponding amplitude (associated with the boundary 
conditions which the scattered field should satisfy) in the two-dimensional impact-parameter- 
space (which is here the x?/-plane), and the integration extends over the region E in which 
a (b) is different from zero. In those cases in which the scatterer is symmetric with respect 
to the scattering axis (here the z-axis), Eq.(|T9"|) can be expressed, by using an integral 
representation for J , as 

oo 

f(q) = — fbdba(b)J (qb) . (20) 

Z7T J 
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where q and b are the magnitudes of q and b respectively. 

The following should be mentioned in connection with Eqs.flTi^D and (PUf): Many of 
the well-known phenomena related to Frauenhofer diffraction have been deducedi! from 
these equations under the additional condition (which is directly related to the boundary 
conditions imposed on the scattered field) \k'\ = \k\ = u' = u, that is, k' differs from k only 
in direction. In other words, the outgoing light wave has exactly the same frequency, and 
exactly the same magnitude of wave-vector as those for the incoming wave. (This means, 
quantum mechanically speaking, the outgoing photons are also on-shell photons, the energies 
of which are the same as the incoming ones.) In such cases, it is possible to envisage that 
q is approximately perpendicular to k and to k', that is, q is approximately perpendicular 
to the chosen z-axis, and thus in the above-mentioned xy-p\ane (that is q ~ q±). While the 
scattering angle distribution in such processes (which are considered as the characteristic 
features of elastic diffractive scattering) plays a significant role in understanding the observed 
diffraction phenomena, it is of considerable importance to note the following: 

(A) Eqs. fll9D and fl20l) can be used to describe diffractive scattering with, or without, this 
additional condition, provided that the difference of k! and k in the longitudinal direction 
(i.e. in the direction of k) is small compared to q± = \q±\ so that q± can be approximated 
by q. In fact, Eqs. (|19"D and fl20|) are strictly valid when <f is a vector in the above-mentioned 
:q/-plane, that is when we write q± instead of q. Now, since Eqs.([T9|) and (p0|) in such 
a form (that is when the replacement q — > q±_ is made) are valid without the condition q 
should approximately be equal to q± and in particular without the additional condition 
\k'\ — \k\ — to' — u, it is clear that they are also valid for inelastic scattering processes. In 
other words, Eqs. ([19]) and ( pD| ) can also be used to describe inelastic diffractive scattering 
(that is, processes in which uj' 7^ u, \k'\ 7^ \k\) provided that the following replacements 
are made. In Eq.©, q q ± , f(q) -> f incl .(q±), a{b) -> a incL (b); and in Eq.(|20|) q -»• q±, 
f(l) ~> /inei.(9J-), <x(b) — > a inel .(6). Hereafter, we shall call Eqs. ([T9|) and (p0|) with these 
replacements Eq. (pjf ) and Eq. (^Of ) respectively. We note, in order to specify the dependence 
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of find, on u' and k'» (that is on w'-w and fen — fey), further information on energy-momentum 
transfer in such scattering processes is needed. This point will be discussed in more detail 
in Section |K]. 

(B) In scattering processes at large momentum-transfer where the magnitude of q± is 
large (|<fj_| 2 0.05 GeV 2 , say), it is less probable to find diffractive scattering events in 
which the additional condition \k'\ = \k\ and uj' = uo can be satisfied. This means, it is 
expected that most of the diffraction-phenomena observed in such processes are associated 
with inelastic diffractive scattering. 

(C) Change in angle but no change in magnitude of wave- vectors or frequencies is likely 
to occur in processes in which neither absorption nor emission of light takes place. Hence, 
it is not difficult to imagine, that the above-mentioned condition can be readily satisfied in 
cases where the scattering systems are time-independent macroscopic apertures or objects. 
But, in this connection, we are also forced to the question: "How large is the chance for 
a incident wave not to change the magnitude of its wave-vector in processes in which the 
scatterers are open, dynamical, complex systems, where energy- and momentum-exchange 
take place at anytime and everywhere?!" 

The picture for inelastic diffractive scattering has two basic ingredients: 
First, having the well-known phenomena associated with Frauenhofer's diffraction and 
the properties of de Broglie's matter waves in mind, the beam particles (7*, 7, p or p 
shown in Fig.|8|) in these scattering processes are considered as high-frequency waves passing 
through a medium. Since, in general, energy- and momentum-transfer take place during the 
passage through the medium, the wave- vector of the outgoing wave differs, in general, from 
the incoming one, not only in direction, but also in magnitude. For the same reason, the 
frequency and the longitudinal component of the wave- vector of the outgoing wave (that is 
the energy, and/or the invariant mass, as well as the longitudinal momentum of the outgoing 
particles) can be different from their incoming counterparts. 

Second, according to the results obtained in Sections [I]-|7| of this paper, the medium is 
a system of color-singlet gluon-clusters (cq's) which are in general partly inside and partly 
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outside the proton — in form of a "cluster cloud" . Since the average binding energy between 
such color-singlet aggregates are of Van der Waals typell, and thus it is negligibly small 
compared with the corresponding binding energy between colored objects, we expect to see 
that, even at relatively small values of momentum-transfer (|£|<1 GeV 2 ,say), the struck Cq 
can unify with (be absorbed by) the beam-particle, and "be carried away" by the latter, 
similar to the process of "knocking out nucleons" from nuclear targets in high-energy hadron- 
nucleus collisions. It should, however, be emphasized that, in contrast to the nucleons in 
nucleus, the Cq's which can exist inside or outside the confinement-region of the proton 
are not hadron-like (See Sections |3|-|6] for more details). They are BTW- avalanches which 
have neither a typical size, nor a typical lifetime, nor a given static structure. Their size- 
and lifetime-distributions obey simple power-laws as consequence of SOC. This means, in 
the diffraction processes discussed here, the size of the scatterer(s), and thus the size of the 
carried- away Cq is in general different in every scattering event. It should also be emphasized 
that these characteristic features of the scatterer are consequences of the basic properties of 
the gluons. 



10. Can such scattering systems be modeled quantitatively? 

To model the proposed picture quantitatively, it is convenient to consider the scattering 
system in the rest frame of the proton target. Here, we choose a right-handed Cartesian 
coordinate with its origin O at the center of the target-proton, and the z-axis in the di- 
rection of the incident beam. The xy-plane in this coordinate system coincides with the 
two-dimensional impact-parameter space mentioned in connection with Eqs.([19f) and fl20f ) 
[which are respectively Eq. (|19D and Eq.(^) after the replacements mentioned in (A) below 



Eq.(p0|)1, while the ?/z-plane is the scattering plane. We note, since we are dealing with in- 
elastic scattering (where the momentum transfer, including its component in the longitudinal 
direction, can be large; in accordance with the uncertainty principle) it is possible to envisage 
that (the c.m.s. of) the incident particle in the beam meets Cq's at one point B = (0,b,z). 
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where the projection of OB along the y-axis characterizes the corresponding impact param- 
eter b. We recall that such Cg's are avalanches initiated by local perturbations (caused by 
local gluon-interactions associated with absorption or emission of one or more gluons; see 
Sections [L] - [7] for details) of SOC states in systems of interacting soft gluons. Since gluons 
carry color, the interactions which lead to the formation of color-singlet gluon-clusters (cq) 
must take place inside the confinement region of the proton. This means, while a consider- 
able part of such Cq's in the cloud can be outside the proton, the location A, where such an 
avalanche is initiated, must be inside the proton. That is, in terms of OA = r, AB = RA{b), 
and proton's radius r p , we have r < r p and [Ra{o)] 2 = b 2 + z 2 + r 2 — 2{b 2 + z 2 ) l ^ 2 r cos LBOA. 
For a given impact parameter b, it is useful to know the distance -R^(fe) between B and A, 
as well as "the average squared distance" (R 2 A {b)) = b 2 + z 2 + a 2 , a 2 = 3/5 r 2 , which is 
obtained by averaging over all allowed locations of A in the confinement region. That is, we 
can model the effect of confinement in cluster-formation by picturing that all the avalanches, 
in particular those which contribute to scattering events characterized by a given b and a 
given z are initiated from an "effective initial point" (Ab), because only the mean distance 
between A and B plays a role. (We note, since we are dealing with a complex system with 
many degrees of freedom, in which B as well as A are randomly chosen points in space, we 
can compare the mean distance between B and A with the mean free path in a gas mixture 
of two kinds of gas molecules — "Species B" and "Species A" say, where those of the latter 
kind are confined inside a subspace called "region p" . For a given mean distance, and a 
given point B, there is in general a set of A's inside the "region p", such that their distance 
to B is equal to the given mean value. Hence it is useful to introduce a representative point 
(A B ), such that the distance between (A B ) and B is equal to the given mean distance.) 
Furthermore, since an avalanche is a dynamical object, it may propagate within its lifetime 
in any one of the Air directions away from (Ab)- (Note: avalanches of the same size may 
have different lifetimes, different structures, as well as different shapes. The location of 
an avalanche in space-time is referred to its center-of-mass.) Having seen how SOC and 
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confinement can be implemented in describing the properties and the dynamics of the Cg's, 
which are nothing else but BTW-avalanches in systems of interacting soft gluons, let us now 
go one step further, and discuss how these results can be used to obtain the amplitudes in 
impact-parameter-space that leads, via Eq.(pO|), to the scattering amplitudes. 

In contrast to the usual cases, where the scatterer in the optical geometrical picture of a 
diffractive scattering process is an aperture, or an object, with a given static structure, the 
scatterer in the proposed picture is an open, dynamical, complex system of Cq's. This implies 
in particular: The object (s), which the beam particle hits, has (have) neither a typical size, 
nor a typical lifetime, nor a given static structure. 

With these in mind, let us now come back to our discussion on the double differential 
cross section d 2 a/dtd(M x l /s). Here, we need to determine the corresponding amplitude 
a incl (6) in Eq.(pOf) [see the discussion in (A) below Eq.fl20|)l. What we wish to do now, is 
to focus our attention on those scattered matter-waves whose de Broglie wavelengths are 
determined by the energy-momentum of the scattered object, whose invariant mass is M x . 
For this purpose, we characterize the corresponding a incl (6) by considering it as a function 
of M x , or Mj/s, or xp. We recall in this connection that, for inelastic diffractive scattering 
processes in hadron-hadron collisions, the quantity M 2 /s is approximately equal to xp, 
which is the momentum fraction carried by the struck Cq's with respect to the incident beam 
(see Fig.[8] for more details; note however that q c , k and p x in Fig.[8] correspond respectively 
to q, k and k! in the discussions here.). Hence, we shall write hereafter a(b\M 2 /s) or 
a{b\xp) instead of the general expression a incl .(&). This, together with Eq. (fffi), leads to 
the corresponding scattering amplitude f(q_i_\xp), and thus to the corresponding double 
differential cross-section d 2 a/dt dxp, in terms of the variables \t\ ~ |<fj_| 2 and xp ~ M x /s in 
the kinematical region: \t\ < M 2 X < s. 
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11. The role played by the space-time properties of the gluon-clusters 

For the determination of a(b\xp), it is of considerable importance to recall the following 
space-time properties of the Cq's which are BTW- avalanches due to SOC: 

(i) SOC dictates, that there are BTW-avalanches of all sizes (which we denote by different 
S values), and that the probability amplitude of finding an avalanche of size 5* can be 
obtained from the size-distribution Ds{S) = S~ IJ ' where the experimental results presented 
in Sections [I] - [7| show: fi « 2. This means, D$(S) contributes a factor S 1-1 , thus a factor 
Xp 1 to the scattering amplitude ot(b \xp). Here, as well as in (ii), we take into account (see 
Sections [| and || for details) that the size S of a Cq is directly proportional to the total 
amount of the energy the Cq carries; the amount of energy is xpP°, where P° is the total 
energy of the proton, and xp is the energy fraction carried by the Cq. 

(ii) QCD impliesil that the interactions between two arbitrarily chosen colored con- 
stituents of Cq is stronger than those between two Cq's, because the latter should be inter- 
actions of Van der Waals type. This means, the struck Cq can unify with the beam-particle 
(maybe by absorbing each other), and viewed from any Lorentz frame in which the beam- 
particle has a larger momentum than the Cq, the latter is "carried away" by the beam 
particle. Geometrically, the chance for the beam-particle to hit an Cq of size S ( on the plane 
perpendicular to the incident axis) is proportional to the area that can be struck by the 
(c.m.s.) of the beam particle. The area is the 2/3-power of the volume S, S 2 ^ 3 , and thus it 
is proportional to x 2 p 3 . 

(iii) Based on the above-mentioned picture in which the Cq's propagate isotropically from 
(Ab), the relative number-densities at different b- values can be readily evaluated. Since for 
a given b, the distance in space between (Ab) and B = (0, b, z) is simply {b 2 + z 2 + a 2 ) 1 / 2 , 
the number of Cq's which pass a unit area on the shell of radius (b 2 + z 2 + a 2 ) 1 / 2 centered at 
{Ab) is proportional to (b 2 + z 2 + a 2 ) -1 , provided that (because of causality) the lifetimes 
(T's) of these Cq's are not shorter than r min (6). The latter is the time interval for a Cq 
to travel from (Ab) to B. This means, because of the space-time properties of such Cq's, 
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it is of considerable importance to note: First, only Cq's having lifetimes T > r min (6) can 
contribute to such a collision event. Second, during the propagation from (Ap) to B, the 
motion of such a Cq has to be considered as Brownian. In fact, the continual, and more 
or less random, impacts received from the neighboring objects on its path leads us to the 
well known0 result that the time elapsed is proportional to the mean-square displacement. 
That is: T mln (b) oc b 2 + z 2 + a 2 . Furthermore, we recall that Cq's are due to SOC, and thus 
the chance for a Cq of lifetime T to exist is Dt{T) oc T~ u where the experimental value (see 
Sections |l| - 0) for v is v ~ 2. Hence, by integrating T~ 2 over T from r min (6) to infinity, 
we obtain the fraction associated with all those whose lifetimes satisfy T > r min (6): This 
fraction is T min (&) -1 , and thus a constant times (b 2 + z 2 + a 2 )~ l . 

The amplitude a(b \xp) can now be obtained from the probability amplitude for Cq- 
creation mentioned in (i), by taking the weighting factors mentioned in (ii) and (iii) into 
account, and by integrating^ over z . The result is: 

a (b \x P ) = const.Xp 1/3 (6 2 + a 2 )~ 3/2 . (21) 

By inserting this probability amplitude in impact-parameter-space, for the beam particle to 
encounter a Cq, which carries a fraction xp of the proton's total energy, in Eq.fl20|') [which is 
Eq.(pOl) with the following replacements: q — > q±, f(q) /meiX^l^p) an d oi{b) ai n oi.(^l x -p)] 
we obtain the corresponding probability amplitude f{q\xp) in momentum-space: 

f(q ± \x P ) = const. f°° bdbxp 1/3 (b 2 + a 2 )~ 3/2 J Q (q±b), (22) 

J 



where q± = \q±\ ~ y\t\ (in the small £p-region, xp < 0.1, say) is the corresponding 
momentum-transfer. The integration can be carried out analyticallyii, and the result is 

f(q±\xp) = const. Xp 1 ^ 3 exp(—aq±) . (23) 

Hence, the corresponding double differential cross-section d 2 a/dtdxp can approximately be 
written as 

-4i- = Nx 'p 2/3 ex p(- 2a #i)> ( 24 ) 

7T dtdxp v 
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where N is an unknown normalization constant. Because of the kinematical relationship 
Xp ~ M 2 /s for single diffractive scattering in proton-proton and proton-antiproton collisions 
(see Fig.||for more details), this can be, and should be, compared with the measured double 
differential cross-sections cPa / dt d(M 2 / s) at different t- and s- values and for different missing 
masses M x in the region M 2 /s <C 1 where q± is approximately J\t\. The comparison is shown 



in Fig.[T(| Here, we made use of the fact that a 2 = 3/5 r 2 , where r p is the proton radius, 
and calculated a by setting r 2 to be the well-knownS mean square proton charge radius, the 
value of which is r 2 = (0.81 fm) 2 . The result we obtained is: a = 3.2 GeV -1 . The unknown 
normalization constant is determined by inserting this calculated value for a in Eq . ([24]) . 
and by comparing the right-hand-side of this equation with the d 2 a/dtd(M 2 /s) data taken 
at \t\ = 0.2 GeV 2 . The value is N = 31.1 mb GeV~ 2 . All the curves shown in Fig.[10| are 
obtained by inserting these values for a and N in Eq.(|24]). 



While the quality of the obtained result, namely the expression given on the right-hand- 
side of Eq.(|2jD together with the above-mentioned values for a and N, can be readily judged 
by comparing it with the data, or by counting the unknown parameters, or both, it seems 
worthwhile to recall the following: The two basic ingredients of the proposed picture which 
have been used to derive this simple analytical expression are: first, the well-known optical 
analogy, and second, the properties of the dynamical scattering system. The latter is what 
we have learned through the data-analysis presented in Sections [l] - [7]. 

Based on the theoretical arguments and experimental indications for the observation (see 
Ref. Jl2 and Section [l] - [7| of this paper for details) that the characteristic features of in- 



elastic diffractive scattering processes are approximately independent of the incident energy 
and independent of the quantum-numbers of the beam-particles, the following results are 
expected: The explicit formula for the double differential cross-section as shown in Eq.(f24"D 
should also be valid for the reactions 7p — > Xp and 7*p — > Xp. While the normalization 
constant iV (which should in particular depend on the geometry of the beam particle) is 
expected to be different for different reactions, everything else - especially the "slope" as 
well as the power of xp should be exactly the same as in pp- and pp-collisions. In this 
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sense, Eg. ([24]) with a 2 = 3/5rp (rp is the proton radius) is our prediction for jp — > Xp and 
7*p — > Xp which can be measured at HERA. 

Furthermore, in order to obtain the integrated differential cross-section da/dt, which has 
also been measured for different reactions at different incident energies, we only need to 
sum/integrate over xp in the given kinematic range (xp < 0.1, say). The result is 

^(t) = Cexp(-2aJ\t\), (25) 

where C is an unknown normalization constant. While this observation has already been 
briefly discussed in the previous LetterEl, we now show the result of a further test of its 



universality: In Fig. 11, we plot 



1 , r 1 da , , . . . 



2^ ° l C dt 

for different reactions at different incident energies in the range 0.2GeV 2 <|t|< 4GeV 2 . Here 
we see in particular that, measurements of da/dt for ^p and jp reactions at larger |£|-values 
would be very useful. 



12. Concluding remarks 

Based on the characteristic properties of the gluons — in particular the local gluon-gluon 
coupling prescribed by the QCD-Lagrangian, the confinement, and the non-conservation of 
gluon-numbers, we suggest that a system of interacting soft gluons should be considered as 
an open dynamical complex system which is in general far away from equilibrium. Taken 
together with the observations made by Bak, Tang and Wiesenfeld (BTW)§0, we are led to 
the conclusion, that self-organized criticality (SOC) and thus BTW-avalanches exist in such 
systems, and that such avalanches manifest themselves in form of color-singlet gluon-clusters 
(cq) in inelastic diffractive scattering processes. 

In order to test this proposal, we performed a systematic data-analysis, the result of 
which is presented in Sections |l] - 0: It is shown that the size-distributions, and the lifetime- 
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distributions, of such cjj's indeed exhibit power-law behaviors which are known as the fin- 
gerprints of SOcEMl. Furthermore, it is found that such exponents are approximately the 
same for different reactions and/or at different incident energies — indicating the expected 
universality and robustness of SOC. Hence, the following picture emerges: For the beam 
particle (which may be a virtual photon, or a real photon, or a proton, or an antiproton; see 
Fig|| for more details) in an inelastic diffractive scattering process off proton (one may wish 
to view this from a "fast moving frame" such as the c.m.s. frame), the target proton appears 
as a cloud of Cg's which exist inside and outside the confinement region of the proton. The 
size-distribution Ds{S) and the lifetime-distribution D?(T) can be expressed as and 
T~ u respectively, where the empirical values for and v are /i^z/^2, independent of the 
incident energy, and independent of the quantum numbers of the beam particles. 

What do we learn from this? Is this knowledge helpful in understanding hadronic struc- 
ture and/or hadronic reactions in Particle Physics? In particular, can this knowledge be 
used to do quantitative calculations — especially those, the results of which could not be 
achieved otherwise? 

In order to demonstrate how the obtained knowledge can be used to relate hadron- 
structure and hadronic reactions in general, and to perform quantitative calculations in 
particular, we discuss the following question — a question which has been with the high- 
energy physics community for many years: 

"Can the measured double differential cross section d 2 a/dtd(M^/s) for inelastic diffrac- 
tive scattering in proton-proton and in antiproton-proton collisions, in the kinematical region 
given by 0.2GeV 2 <|t|< 3.25 GeV 2 , 16GeV<v^< 630 GeV, and M 2 /s <0.1, be understood 
in terms of optical geometrical concepts?" 

The answer to this question is "Yes!", and the details are presented in Sections |9| 
- [11] where the following is explicitly shown: The characteristic features of the existing 
d 2 a/dt d(M 2 /s)-data are very much the same as those in optical diffraction, provided that 
the high-energy beams are considered as high-frequency waves, and the scatterer is a system 
of color-singlet gluon-clusters (cq) described in Sections [l] - [7| of this paper. Further mea- 
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surements of double differential cross sections, especially in 7*p- and 7p-reactions, will be 
helpful in testing the ideas presented here. 
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k k' 




Fig. 1b 

Fig. 1. The well-known Feynman diagrams (a) for diffractive and (b) for normal deep-inelastic 
electron-proton scattering are shown together with the relevant kinematical variables which describe 
such processes. 
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Fig. 2. F^ {3) (P,Q 2 ;x P )/P is plotted as a function of j3 for given xp-intervals and for fixed 
Q 2 -values. The data are taken from Ref.||. The lines are only to guide the eye. 
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Fig. 3. (a) F 2 D(3) (/3,Q 2 ;x P )/x P is plotted function of xp for different values of /3 and Q . 
The data are taken from Ref.||. (b) same as (a) but with data taken from Ref.[|j. Note that 
in these log-log plots, almost all existing data points lie on straight-lines with approximately the 
same slope irrespective of the values of Q 2 and/or (3. 
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Fig. 4. F2 {3 \f3,Q 2 ;x P )/f3 is plotted function of xb in the indicated (5- and Q -ranges. 



The data are taken from Ref.Q. 
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Fig. 5b 



Fig. 5. (a) F 2 D(3) (/?,Q 2 ;x B //3)/(/?x B ) is plotted function of xb for fixed (3- and Q 2 -values. 
The data are taken from Ref.||. (b) same as (a) but with data taken from Ref.[|j. Note that 
in these log-log plots, almost all existing data points lie on straight-lines with approximately the 
same slope irrespective of the values of Q 2 and/or j3. 
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(1/Xp)d 2 c D /dQ 2 dx p (0.389 mb/GeV 2 ) 
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slope: -2.0210.23 



Fig. 6. (a) (l/xp)d 3 a D /df3dQ 2 dx P is plotted as a function of xp in different bins of /3 and 
Q 2 . The data are taken from Ref.g. (The factor 0.389mb is due to (lGeV) -2 = 0.389mb.) (b) 
(l/xp)d 2 a D /dQ 2 dxp in units of GeV -4 is plotted as a function of xp in different bins of Q 2 . The 
data are taken from Ref.Q. (c) same as (a) but with data taken from Ref.Q]. (The factor 0.389 mb 
is due to (lGeV)- 2 = 0.389 mb.) 
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(1/Xp)dc/dx p (jib) 



10 10 
10 9 


• 


H1 data 




T D\ 

r ^ \ • 


. <W>=187GeV, M Y <1.6GeV 
. <W>=187GeV, M Y >1.6GeV 
o <W> =231 GeV, M Y < 1 .6 GeV 
□ <W> =231 GeV, M Y > 1 .6 GeV 


10 8 




ZEUS data (normalized to c tol ) 
a <W>=200GeV i 


10 7 


r 






10 6 


r 






10 5 








10 4 








10 3 








10 2 


r power=1.98±0.07 






-5 -4 
10 10 


10 " 3 


10 2 



Fig. 7 



x p =l\/|2/W 2 



Fig. 7. {1 / x p)da / dx p for photoproduction ^ + p ^ X + p \s plotted as a function of xp. The 
data are taken from Ref.[^]. Note that the data in the second paper are given in terms of relative 
cross sections. Note also that the slopes of the straight-lines are the same. The two dashed lines 
indicate the lower and the upper limits of the results obtained by multiplying the lower solid line 



by a^Qi = 154 ± 16(stat.) ± 32(syst.)//6. This value is taken from the third paper in Ref.[21]. 
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Fig. 8. Diagrams for different single diffractive reactions, together with the definitions of the 
relevant kinematic variables. 
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Fig. 9. a) (1/ xp)d 2 a / dxpdt for single diffractive p + p — ► p + X and p + p —> p + X reactions 
is plotted as a function of xp at different values of t and y/s. The data are taken from Ref.|5[[|. 
b) The integrated (with respect to two different |t|-ranges) differential cross section (l/xp)da/dxp 
for single diffractive p + p — > p + X and p + p ^ p + X reactions is plotted as a function of xp. 
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Fig. 10. The double differential cross section (1/ir) dPa/dt d(M% / s) for single diffractive pp 
and pp reactions is shown as function of xp at fixed values of t where 0.15 GeV 2 <|t|<3.25 GeV 2 . 
The data are taken from Refs. ||^]. The solid curve is the result obtained from Eq.(p4|). The 
dashed curve stands for the result obtained from the same formula by using the t-value given in 
the bracket. 
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Fig. 11. The quantity (— l/(2- v /j*[)) log da/dt] is plotted versus \/\t\ for different single 
diffractive reactions in the range 0.2GeV 2 <|t|<4GeV 2 . The data are taken from Refs. igUgl]. 
Here, C, the normalization constant is first determined by performing a two-parameter fit of the 



corresponding da/dt-d&t& to Eq.(25). 
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